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INTRODUCTION 

A great  deal  of  current  interest  in  piezoelectric  vibrators  focuses  at 
the  low  end  of  the  frequency  spectrum,  primarily  for  timing  applications. 

In  the  solution  of  the  vibration  problem  for  the  bars  and  forks  used  at  low 
frequencies,  it  is  usually  convenient  to  employ  the  elastic  compliance 
matrix  s instead  of  the  stiffness  matrix  c.’>^ 

In  those  situations  where  the  stiffnesses  c are  used,  e.g. , the  high 
frequency  thickness  modes  of  plates,  the  inclusion  of  acoustic  loss  is 
accomplished  by  allowing  c to  become  complex.^  The  material  coefficients 
appearing  in  the  imaginary  part  are  the  elements  of  the  viscosity  matrix  n- 
In  the  sequel  we  give  the  corresponding  fluency  matrix  arising  from  the 
complex  compliance,  and  compute  the  elements  for  quartz. 

FLUENCY  MATRIX 

The  influence  of  internal  friction  upon  the  elastic  constants  c and  s 
was  discussed  briefly  by  Voigt. '*  Under  suitable  conditions  of  approximation 
Voigt  gives,  in  our  notation,  the  Hooke's  law  relations  in  the  presence  of 
loss; 

T = c*S  = (c  + n9/3t)S,  (1) 

S = S*T  = (s  - t9/9t)T.  (2) 

The  compliance  matrix  may  be  developed  from  these  equations  as  a power  series 
in  c and  q.  For  harmonic  time  dependence  one  obtains 

2n~ 

(-)"  0)2"  (nc-l) 

2n  — I 

(-)"  0)2"  (c-lq)  jc-\  (3) 

and 

C = c‘^ns.  (4) 


1.  W.  G.  Cady,  Piezoelectricity.  New  York:  McGraw-Hi 1 1,  1946,  and  Dover, 1964. 

2.  "Standards  on  Piezoelectric  Crystals,  1949."  Proc.  IRE,  Vol . 37, 
pp.  1378-1395,  Dec.  1949.  (IEEE  Standard  No.  176.) 

3.  J.  Lamb  and  J.  Richtjr,  "Anisotmnic  Acoustic  Attenuation  with  New 
Measurements  for  Quartz  at  Room  lemperature,"  Proc.  Roy.  Soc.  (London), 
Vol.  293A,  pp.  479-492,  1966. 

4.  W.  Voigt,  Lehrbuch  der  Kristallphysik.  Leipzig:  Teubner,  2nd  ed.,  1928, 

pp.  792-796. 
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The  analogous  inverse  relations  are  obvious. 


Alternatively,  one  may  write 

c*s*  = (c  + j(jjn)(s  - jwrj  = I,  (5) 

multiply  out,  equate  real  and  imaginary  parts,  and  solve  to  get: 

s = (c  + w^nc'^n)  , (6) 

A = c'^ns,  (7) 

with  corresponding  inverse  relations.  In  the  limit  of  zero  frequency. 
Equations  (3)  and  (6)  become 

s = c'^ , (8) 

and  Equations  (4)  and  (7)  are 

C = c“^nc”^  = sns.  (9) 

APPLICATION  TO  QUARTZ 

Using  the  room  temperature  values  for  c^  from  Bechmann^  with  the  n 
values  of  Lamb  and  Richter^  in  Equation  9 gives  the  values  for  the  elements 
of  the  fluency  matrix  as  shown  in  Table  1.  The  c matrix  has  the  same  symmetry 
as  s for  all  crystal  classes.  Recalculating  ^ with  w = lO^^  changes  the 
entries  by  a few  units  in  the  third  decimal  place.  Table  2 contains  the 
fluency  matrix  elements  for  w = 10^^. 

TABLE  1. 

FLUENCY  MATRIX  FOR  QUARTZ  AT  w = 0 


1 1 

18.56 

Cl,  4 

17.65 

3 3 

6.21 

Cg6 

30.70 

1 7 

3.21 

Cl  4 

7.87 

1 3 

3.80 

(all 

in  10-^6  s/Pa 

5.  R.  Bechmann,  "Elastic  and  Piezoelectric  Constants  of  Alpha-Quartz," 
Phys.  Rev.,  Vol . 110,  pp.  1060-1061,  June  1958. 
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TABLE  2 


FLUENCY  MATRIX  FOR  QUARTZ  AT  id  = 10^^ 


5.30 

9.01 

1 1 

4 4 

3.42 

13.52 

3 3 

66 

-1.46 

c 

2.56 

1 ? 

1 4 

-0.46 

1 3 

(all 

in  10“^®  s/Pa) 

MOTIONAL  TIME  CONSTANTS 


From  the  solution  to  a particular  vibration  problem,  the  eigenvalue 
or  C[^*  for  the  mode  m in  question  will  arise.  Motional  time  constants 

t('^)  and  rCfti)  are  then  defined  as 
1 1 


T 


(m) 


Im(c^*)  rin, 
Lim  = 

0)^  o)Re(c^*) 


(10) 


and 


. . ^m 

Lim 

ur*-o  o)Re(s  *)  s„ 

' m m 


(11) 


Application  of  Equation  (10)  to  doubly  rotated  quartz  plates  has  recently 
been  made,®  where  it  is  shown  how  may  be  obtained  from  n and  from  the 
solution  to  the  lossless  problem. 

xjm)  is  equivalently  expressed  in  terms  of  the  constants  of  the 
equivalent  circuit  as 

= R C = 1 (12) 

1 1 1 TTrr  ‘ 


Additional  relations  may  be  given  in  terms  of  acoustic  attenuation  or  of 
logarithmic  decrement. 

Equation  (12)  permits  calculation  of  the  intrinsic  Q of  a mode, 
exclusive  of  extrinsic  factors  such  as  mounting  loss  and  ambient  loading. 
It  also  pertains  to  f(l’)  for  bar  and  fork  resonators. 


6.  A Ballato,  "Doubly  Rotated  Thickness  Mode  Plate  Vibrators,"  in  Physical 
Acoustics:  Principles  and  Methods  (W.  P.  Mason  and  R.  N.  Thurston,  eds.) 

New  York:  Academic,  1977,  Vol.  13,  pp.  115-181. 
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For  quartz  at  room  temperature,  t and  t run  about  10  fs,  with  corres- 
ponding relaxation  frequencies  falling'in  the  I’nfrared  region.  It  is  no 
accident,  therefore,  that  the  IR  absorption  spectrum  is  used  as  a quality 
control  criterion  for  inspection  of  cultured  quartz’* the  absorption  is 
found  to  be  correlated  with  the  intrinsic  Q of  resonators  made  from  the 
material . 

Very  recently,  a new  procedure  has  been  proposed  for  evaluating  the 
quality  of  cultured  quartz  samples.’®  The  method  requires  resonator  Q 
measurements  made  on  extensional  mode  quartz  bars.  The  data  lead  to  experi- 
mental values  of  Tj  for  comparison  with  Equation  (11).  The  ratio  of  measured 
Tj  to  calculated  furnishes  the  condition  of  acceptability. 

A constraint  on  the  method  is  the  Q-limiting  effect  of  the  mounting 
supports  of  the  resonator  bar.  Sherman’ ‘ has  given  an  approximate  relation 
that  prescribes  limits  on  the  dimensions  of  the  mounting  supports  and  on  the 
accuracy  of  their  placement  about  the  bar  nodal  positions  to  ensure  the 
dominance  of  the  intrinsic  Q. 

THE  -18.5°  QUARTZ  BAR 

Mason  introduced  the  -18.5°  quartz  bar.”  Its  orientation  is  denoted’ 
as  (XYt)ij;  = -18.5°,  and  is  shown  in  Figure  1.  This  cut  vibrates  in  an  exten- 
tional  mode  along  its  length.  Its  angle  of  cut  is  determined  by  location  of 
the  zero  of  the  elastic  compliance  s'  : 

? 4 


7.  D.  M.  Dodd  and  D.  B.  Fraser,  "The  3000-3900  cm"^  Absorption  Bands  and 
Anelasticity  in  Crystalline  a-Quartz,"  J.  Phys.  Chem.  Solids,  Vol . 26, 
pp.  673-686,  1965. 

8.  D.  B.  Fraser,  D.  M.  Dodd,  D.  W.  Rudd,  and  W.  J.  Carroll,  "Using  Infrared 
to  Find  the  Mechanical  Q of  a-Quartz,"  Frequency,  Vol.  4,  pp.  18-21, 

Jan. -Feb.  1966. 

9.  J.  Asahara  and  S.  Taki,  "Physical  Properties  of  Synthetic  Quartz  and  Its 
Electrical  Characteristics,"  in  Proc.  26th  Annu.  Frequency  Control  Symp., 
US  Army  Electronics  Command,  Fort  Monmouth,  NJ,  June  1972,  pp.  93-103. 

10.  H.  Fukuyo,  N.  Oura,  and  F.  Shishido,  "A  New  Quality  Evaluation  Method  of 
Raw  Quartz  by  Measuring  the  Q-Value  of  Y-Bar  Resonator,"  in  Proc.  31st 
Annu.  Frequency  Control  Symp.,  US  Army  Electronics  Command,  Fort  Monmouth, 
NJ,  June  1977,  in  press. 

11.  J.  H.  Sherman,  Jr.,  private  communication,  Aug.  1977. 

12.  W.  P.  Mason,  "Electrical  Wave  Filters  Employing  Quartz  Crystals  as 
Elements,"  Bell  Syst.  Tech.  J.,  Vol.  13,  pp.  405-452,  1934. 
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Figure  1.  Orientations  of  AT-plate  and  -18.50  Quartz  Bar 
s^,  = {-2s  cos^tjj  sini|;  + 2s  sin^ij; 

+ (s^_^  + 2Sj^)  sintp  (cos^ijj  - sin^ij;) 

+ Sj^  cosi|j  (3  sin^ijj  - cos^i|;)}*  cosij;. 

The  compliance  governing  the  frequency  is  s^^: 

s^  = s cos‘*ip  + s sin'*!];  + (s  + 2s  ) cos^tl)  sin^ii) 

-2s  cos^iJi  sinOi. 

14 

The  bar  is  driven  by  an  x-directed  electric  field;  the  piezoelectric 
constant  is: 

/ 

dj2  = -djj  cos^ij;  + dj^  cosi|/  siniji, 
and  the  dielectric  permittivity  is 


(15) 
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Although  the  vibrator  is  referred  to  as  the  "-18.5°  quartz  bar,"  the  more 
recent  elastic  values  of  Bechmann®  show  that  s'^^  = 0 at  = -19.70°.  At 
this  angle,  the  appropriate  numerical  constants  describing  the  resonator  are: 

• s'  = 14.46  X 10-^2  Pa'^ 

2 2 

• N = (4ps'”  = 2.555  MHz-mm 

2 2 

• d'  = -2.278  pC/N 

12 

• £^2  = 39.97  pF/m 

• r = C /C  = TrV(8k2)  = 137.4 

0 1 

• = ^22/^  = 290.9  fF/m. 

Substitution  of  for  sij  in  Equation  (14)  converts  it  into  the 
formula  for  The  entries  from  Table  1 are  then  used  to  compute  the 

effective  value: 

• = 21.63  10-26  s/pa. 

This  leads  to  the  additional  constants: 

• = ^2^/<2  = 

• Pj  = = 5.14  X 10"2  n-m. 

A relaxation  time  constant  appears  from  other  considerations  in  a letter 
by  Cook  and  Wasilik'^  on  the  -18.5°  quartz  bar. 

In  order  to  put  the  size  of  into  perspective,  it  is  instructive  to 
observe  that  in  the  time  Tj  the  longitudinal  elastic  waves  traveling  along 
the  bar  length  coveir  a distance  of  0.76  S,  or  only  a little  in  excess  of  a 
Bohr  radius  (0.529  A).  The  remaining  quantities  appearing  above  are  related 
to  the  parameters  of  the  lumped  electrical  equivalent  circuit,  and  are 
defined  in  Bechmann's  paper. 

Figure  2 shows  the  analog  electrical  equivalent  circuit  superimposed  on 
a sketch  of  the  bar.  Piezoelectric  excitation  takes  place,  according  to 


13.  R.  K.  Cook  and  J.  H.  Wasilik,  "Anelasticity  and  Dielectric  Loss  of 
Quartz,"  J.  Appl . Phys.,  Vol  27,  pp.  836-837,  1956. 

14.  R.  Bechmann,  "Schwingkristal le  fiir  Siebschaltungen,"  Arch.  d.  Elektr. 
Ubertrag.,  Vol.  18,  pp.  129-136,  Feb.  1964. 


k 


6 


\ / 


Figure  2.  Analog  Equivalent  Circuit  of  Extensional  Mode  Bar. 


this  scheme,  at  the  bar  ends.  The  transformers  are  accordingly  placed  there 
The  input  admittance  is'^ 


Yin  = Zin  = > ^0  ^ ^ ’ 


15.  D.  A.  Berlincourt,  D.  R.  Curran,  and  H.  Jaffe,  "Piezoelectric  and 

Piezomagnetic  Materials  and  Their  Function  in  Transducers,"  in  Physical 
Acoustics:  Principles  and  Methods  (W.  P.  Mason,  ed.)  New  York: 


Academic,  1964,  Vol . TA,  pp.  169-270. 
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CONCLUSIONS 

We  have  shown  how  the  acoustic  loss  at  high  frequency  may  be  used  to 
derive  an  associated  set  of  constants  for  low  frequency  use.  These  have  been 
denoted  as  the  fluencies,  have  been  calculated  for  quartz,  and  have  been 
applied  to  the  vibrations  of  bars.  The  low  frequency  time  constant  Tj  has 
been  introduced,  as  well  as  the  analog  equivalent  circuit  for  low  frequency 
resonators.  A useful  future  area  of  work  consists  in  measuring  the  tempera- 
ture coefficients  of  the  viscosities  and/or  fluencies. 
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